NBI-HE-00-45 



1 



o 
o 
o 

> 

o 

(N 



> 

o 



o 
o 

in 



X 



Non-equilibrium field theory* 

Dietrich B6deker a t 

a The Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 

I discuss various topics in relativistic non-equilibrium field theory related to high energy physics and cosmology. 
I focus on non-perturbative problems and how they can be treated on the lattice. 



1. INTRODUCTION 

The subject of non-equilibrium field theory is 
the evolution of many-particle systems in real 
(Minkowski) time. It also gives new insights into 
the much better understood equilibrium proper- 
ties of quantum fields. Most of its applications in 
high energy physics are related to the early uni- 
verse and to heavy ion collisions. 

Many interesting problems in non-equilibrium 
field theory are non-perturbative. One could hope 
that they can be put on a lattice and then be 
solved by a computer. However, lattice simula- 
tions of quantum field theory (QFT) work in Eu- 
clidean time and nobody knows how to perform 
real time simulations. 

A considerable simplification occurs for sys- 
tems which are in partial (or incomplete) equilib- 
rium (see, e.g., 0). It means that some slowly 
relaxing quantities y a differ significantly from 
their equilibrium values, while all other degrees 
of freedom are thermalized. The most extreme 
example are hydrodynamic modes, their relax- 
ation times diverge when their wavelengths go 
to infinity (other examples will be discussed in 
Sec. ||) . The non-equilibrium problem can then be 
reduced to evolution equations for the y a . These 
equations contain so called transport coefficients, 
and the task of non-equilibrium field theory is to 
determine them from the underlying QFT. 

Partial equilibrium is a good approximation for 
most of the history of the early universe (an im- 
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portant exception is (p)reheating at the end of in- 
flation, that is the transition to a radiation dom- 
inated epoch, see Sec. 2d). Needless to say that 
the most solid results have been obtained for this 
case. 

Fluctuation dissipation relations allow to ex- 
press transport coefficients in terms of thermal 
averages like 



([O(i),O(0)]>, 



(1) 



where the angular brackets denote the average 
over a thermal ensemble with temperature T, 
-H/T\. 



(...)= Z~hr(. 



(2) 



H is the Hamiltonian. For example, electric con- 
ductivity is obtained from (|l|) with O being the 
electromagnetic current (color conductivity, on 
the other hand, has no interpretation in terms 
of a gauge invariant correlation function of some 
current. Instead it arises as a "Wilson coefficient" 
in an effective theory for long distance modes of 
non-abelian gauge fields, see Sec. 4.2). 

I have already mentioned that a serious diffi- 
culty of thermal field theory is the fact that one 
has to deal with real time. Let me illustrate this 
for the case of the expectation value ([]]). One 
could think of expanding Eq. (|lj) in powers of t. 
Then one can compute every coefficient in the 
expansion from the usual Euclidean path integral 
of thermal QFT. But transport coefficients are 
determined by the small frequency limit of the 
Fourier transform of (]]]). This depends on the 
behaviour for t — > oo which can not be captured 
by an expansion around t = 0. 

Fortunately there is a limit of quantum field 
theory in which a non-perturbative treatment is 
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possible. It is the classical field limit of bosonic 
fields. It applies when the number of quanta in 
the field modes of interest is large. This is in- 
deed the case in a variety of interesting prob- 
lems. I discuss the classical field approximation, 
together with some recent applications in Sec. |^. 
A wide class of approximations, many of them 
related to a large N expansion, which are not re- 
stricted to the classical field limit is discussed in 
Sec. [|. In Sec. [| I consider systems in incomplete 
equilibrium. I show how one can use the classical 
field approximation even when quantum effects 
are important by constructing effective classical 
field theories. Sec. || is a brief summary of my 
talk. 

2. CLASSICAL FIELD APPROXIMA- 
TION 

While quantum systems are impossible to sim- 
ulate in real time, there is no principle obstacle to 
performing real time simulations in classical field 
theory. All one has to do is to solve classical field 
equations of motion with the appropriate initial 
conditions, and determine the physical quantity 
of interest from the solution. 

A well known special case of the classical field 
approximation (usually not referred to as such) is 
the dimensional reduction of a d + 1 dimensional 
thermal QFT in the imaginary time formalism to 
the d dimensional theory for the zero Matsubara 
frequency modes of the bosonic fields Q]. This 
is the formal classical limit because for T% — > 
the non-zero modes become infinitely heavy and 
decouple. The thermodynamics of QFT can be 
studied with 4-dimcnsional lattice simulations, 
which is the only reliable way to access the prop- 
erties of hot QCD near the critical temperature. 
Dimensional reduction is convenient for very high 
temperatures when there is a large separation of 
the scale 2irT and the screening length(s) of the 
system. In non-equilibrium field theory the role 
of the classical field approximation is much more 
important since there is no tractable analogue of 
the 4-dimensional Euclidean theory. 

The classical field approximation should be re- 
liable when the number of field quanta in each 
relevant field mode is large. There are indeed in- 



teresting applications where this is the case. The 
remainder of this section is a list of some of them. 

2.1. Preheating after inflation 

Inflation is the only known solution to the hori- 
zon and flatness problem in standard Big Bang 
cosmology. Furthermore, it generates density 
fluctuations which can seed the structure forma- 
tion once electrons and nuclei have combined to 
form atoms. There has been a lot of interest in 
the transition from an inflationary to the radia- 
tion dominated epoch after it was realized || that 
the initially homogeneous inflaton field if can de- 
cay very rapidly into low momentum modes of the 
inflaton itself or into modes of other scalar fields x 
through parametric resonance. This mechanism 
is referred to as pre-heating, and it has drastically 
changed the picture of reheating after inflation. 

The large amplitudes of the amplified modes 
make this problem non-perturbative. At the same 
time, the occupation numbers of these modes 
grow very large which opens the possibility to 
study this problem in the classical field approx- 
imation on the lattice |4|J^]. As an illustration, 
Fig. 1 shows the occupation number of the x field 
(taken from Ref. §). 

Lattice simulations have shown that non-linear 
effects play an important role in reheating [^|J^]. 
Effects like non-thermal phase transitions and de- 
fect formation have been investigated [@-|l(| . Pre- 
heating and the possibility of non-thermal phase 
transitions in gauge theories were studied in . 

2.2. Small x gluon distributions 

The gluon density in hadrons and nuclei be- 
comes large at small x, where x is the fraction 
of the momentum of the hadron/nucleus which 
is carried by the gluon. Then the usual per- 
turbative treatment of parton evolution breaks 
down. In jl2| a model was proposed which de- 
scribes the small x gluon field of large nuclei clas- 
sically. It was first applied to nuclear collisions 
in Ref. jl3). If one is interested in the produc- 
tion of gluons with small transverse momenta a 
non-perturbative treatment of this model is nec- 
essary. A lattice version was developed in Q. 
It was used to solve the non-linear equations of 
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Figure 1. Occupation numbers of x-modes vs 
momentum. The curves correspond to differ- 
ent values of time. There are three visible res- 
onances caused by the oscillating inflaton back- 
ground (from Ref. [6]). 



motion for the gluon fields of two colliding nuclei, 
and energy densities and gluon multiplicities in 
the final state were extracted fl5[| . 

2.3. Tests of approximation schemes 

The approximation schemes to be discussed in 
Sec. [| are usually formulated in quantum field 
theory. It is difficult to asses their reliability. In 
Ref. |L6|] it was pointed out that they can equally 
well be formulated in the classical field limit. It is 
then possible to perform lattice simulations and 
compare the "exact" lattice results with the ap- 
proximate ones. For a scalar model in (1+1) di- 
mensions it was found that early-time behaviour 
is reproduced qualitatively by the Hartree ap- 
proximation, and that including scattering im- 
proves this to the quantitative level. At late times 
the lattice system thermalizes, but the improved 
Hartree approximation fails to reproduce this be- 
haviour |l7| . 

2.4. Droplet nucleation in a thermal first 
order phase transition 

Starting from a high temperature, a first order 
phase transition proceeds through the nucleation 
of droplets, which then grow and eventually fill 



the whole space with the low temperature phase. 
The droplet nucleation is caused by thermal fluc- 
tuations. Most calculations of the nucleation rate 
use the formalism due to Langer [Q which is 
based on a saddle point approximation around 
the the critical droplet. A lot of lattice compu- 
tations have been performed, most of them have 
considered a moderately strong first order phase 
transition (for recent results, see |l9|]). Ref. |20| 
has considered a very strongly first order tran- 
sition in the electroweak theory. The results 
were compared with various analytic calculations 
based on Langer's approach. It was found that 
the accuracy of the analytic results strongly de- 
pend on the approximation used for the critical 
droplet action. The best agreement was found 
when the 2-loop effective potential for the Higgs 
field |2l];|2^] together with the field dependent 
wave function renormalization |2^,|24| is used. 

3. BEYOND THE CLASSICAL FIELD 
APPROXIMATION 

The classical field approximation is clearly re- 
strictive, and it is not capable of treating genuine 
quantum effects. For example, the thermaliza- 
tion of the universe after preheating (cf. Sec. |^) 
is not classical. It is sometimes possible to obtain 
effective classical theories when the non-classical 
modes are weakly interacting and can be inte- 
grated out in perturbation theory (see Sec. 4.2). 
In general, however, other approximations are 
necessary. 

The Hartree approximation was discussed in 
the talk by J. Vink [^5| and there are several re- 
lated approaches. Usually the basic idea is to 
work directly with non-equilibrium Green func- 
tions or, alternatively, with the generating func- 
tional of 1-PI Green functions (g6[ . These satisfy 
certain Schwinger-Dyson-like equations which are 
an infinite set of equations containing the full 
set of n-point functions. To make this prob- 
lem tractable, the hierarchy of Schwinger-Dyson 
equations must be truncated. It is possible to 
formulate such truncations as a systematic expan- 
sion in 1 /N, where N is the number of some fields 

mi- 

Recent lattice studies have focused on the ques- 
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tion whether one can obtain thermalization star- 
ing from some non-equilibrium initial state (see 
also Sec. |2.3| ). Thermalization can be obtained in 
the Hartree approximation if the mean field is al- 
lowed to be inhomogeneous, which was discussed 
in J. Vink's talk |^5| . Thermalization of fermionic 
modes, which can never be described by classical 
fields, was also observed in the inhomogeneous 
mean field approximation p9| . 

Rcf. considers (1+1) dimensional A</? 4 the- 
ory in the loop expansion of the generating func- 
tional of 2-particle irreducible non-equilibrium 
Green functions. At 3-loop order, when scatter- 
ing is taken into account, thermalization was ob- 
served. 

4. CLOSE TO THERMAL EQUILIB- 
RIUM 

As I mentioned in the introduction, there are 
interesting non-equilibrium processes in which 
most of the degrees of freedom are equilibrated 
and only some slowly relaxing ones are not. The 
task of non-equilibrium field theory is then to de- 
termine the relevant diffusion coefficients. 

In weakly coupled theories it is possible to com- 
pute some transport coefficients using perturba- 
tion theory. This is not an easy task. Typi- 
cally one has to sum an infinite set of diagrams, 
which was done for the shear and bulk viscos- 
ity in (3+1) dimensional scalar theory with cu- 
bic and quartic self interactions (3^] . The result 
of |5l| is equivalent to the solution of the clas- 
sical Boltzmann equation with appropriate coef- 
ficients. Other transport coefficients were calcu- 
lated starting directly from the Boltzmann equa- 
tion, which is technically much simpler than a di- 
agrammatic analysis |32| . Another application of 
the Boltzmann equation of current interest is the 
generation of a lepton asymmetry due to heavy 
Majorana neutrinos 

Unfortunately, it is not clear whether the clas- 
sical Boltzmann equation is the leading order ef- 
fective theory in some systematic weak coupling 
expansion which would allow to go beyond the 
leading order. There are problems where this is 
necessary, like for example in electroweak baryo- 
genesis Q . 



There are transport coefficients which are not 
computable in (appropriately resummed) pertur- 
bation theory even when the coupling constant 
is small. A famous example is the rate Jb+l 
for baryon (B) plus lepton number (V) dissipa- 
tion in the standard electroweak theory J34| . The 
Lagrangian of the electroweak theory has a U(l) 
symmetry which would imply that B + L is con- 
served. However, this symmetry is anomalous, 
and B + L violating processes are possible. Sup- 
pose there is some non-zero density Nb+l in the 
universe. It then gets washed out according to 
d 

dj. N B+L = -1b+lNb+l- (3) 

There is a fluctuation dissipation relation stat- 
ing that 7b+l is proportional to the Chern- 
Simons diffusion rate Yes of the weak SU(2) 
gauge fields. Tcs is dominated by soft gauge 
field modes with typical momenta of order g 2 T . 
This is the so called magnetic screening scale at 
which perturbation theory for hot Yang-Mills the- 
ory breaks down ]3q] . Fortunately the smallness 
of this scale allows the use of the classical field 
approximation to compute Teg- The occupation 
number of the bosonic field modes can be esti- 
mated from the Bose-Einstein distribution func- 
tion 

"( k )=dW^T"M (l k ' <<T ) ^ 

which is of order 1/g 2 for the soft (|k| ~ g 2 T) 
modes. 

The question how to calculate Tqs has a long 
history. One reason for this long lasting interest 
is that the B+L dissipation rate is a fundamental 
quantity of the Standard Model of electroweak in- 
teractions (SM) and extensions thereof. It plays 
an important role in particle physics scenarios for 
generating the baryon asymmetry of the universe. 
The use of classical lattice field theory for calcu- 
lating r C s was first suggested in Ref . j$6) , where 
a (1+1) dimensional model was studied. 

Originally it was assumed that the naive clas- 
sical field approximation correctly describes the 
dynamics of the soft modes in (3+1) dimensional 
Yang-Mills theory. Therefore the lattice version 
of the classical Yang-Mills equations of motion 

£> E + DxB = 0, (5. a) 



D E 



(5.b) 



were solved and used to determine the Chcrn- 



Simons diffusion |3q-|40 
E and B are the color 



(see also [J47)). Here 
electric and magnetic 



fields and denotes the covariant derivative 
+ gA^ in the adjoint representation. The 
above assumption would be correct if high mo- 
mentum (|k| ^> g 2 T) modes decouple from the 
soft dynamics. It is suggestive because it is true 
for the thermodynamics of hot Yang-Mills fields: 
At leading order in g the soft non-perturbative 
modes are described by the 3-dimensional pure 
Yang-Mills theory of dimensional reduction, and 
the only role of high momentum modes is to make 
the coupling constant equal to its renormalizcd 
value. 

4.1. Problems of the classical field approx- 
imation 

Classical thermal field theory as a continuum 
theory does not exist. ^ This has been known 
for more than 100 years and it has played a cru- 
cial role in the discovery of quantum theory by 
Planck. 

In dimensional reduction the UV divergences 
are taken care of by adding the appropriate local 
counterterms to the 3-dimensional action. Then 
long distance (I 3> T -1 ) correlation functions of 
the 4-dimensional theory are correctly reproduced 
by the 3-dimensional one. 

However, it was realized in [[ll| that there are 
UV divergences in classical thermal Yang-Mills 
theory which do not occur in equal time corre- 
lation functions []. These divergences are closely 
related to what is called "hard thermal loops" in 
thermal QFT |52] |. They are particularly ugly, 
they are non-local in space and time and on the 
lattice they are sensitive to the lattice geome- 
try (see also Q). The occurrence of such di- 
vergences in lattice field theory is possible due 



3 I use "color" as a generic term for some non-abelian gauge 
charge, not as something specific to QCD. 
4 Take for example the energy density of a free massless 
bosonic field at temperature T. It is proportional to TA 3 , 
where A is the UV cutoff. 

5 Real time correlation functions in scalar Ai/? 4 theory, on 
the other hand, are rcriormalized by the counterterms of 
dimensional reduction 1421 . 
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Figure 2. One- loop polarization tensor in non- 
abelian gauge theory. Only the first two diagrams 
contribute to the UV divergence in the classical 
field limit, while the ghost loop contribution is 
suppressed. 



to the lack of Lorentz invariance which is explic- 
itly broken by the presence of the thermal bath. 
As a consequence, Green functions with external 
momenta p a depend not only on the Lorentz in- 
variants p a ■ pp but on frequencies and spatial 
momenta separately. 

Since such divergences are probably not famil- 
iar to most of the audience, let me describe them 
in some more detail. They are due to the one-loop 
polarization tensor (see Fig. 2). First consider the 
case that the external momentum k is purely spa- 
tial, fco — 0. This is what one encounters in the 
3-dimensional reduced theory. The result is of the 
form 



n^(0,k) - g 2 T 



A d 3 q fk-q k 2 



(6) 



\ q- 

where the function / is dimensionless, and A is 
the UV cutoff. / can be expanded in powers of its 
arguments, and only the lowest order term /(0, 0) 
gives rise to a linearly divergent integral, 



IP^O, k) - g 2 TA/(0, 0) + UV finite. 



(7) 



Gauge invariance ensures that this divergence oc- 
curs only in LI 00 . This is nothing but the familiar 
divergence of the mass term for the adjoint Higgs 
field Aq which is taken care of by adding a mass 
counterterm to the 3-dimcnsional action. 

Now consider the analogous contribution when 
the external frequency is non-zero. In the clas- 
sical field limit the result for IF" is similar to 
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Eq. (|6|). However, the function / now depends on 
three dimensionless variables, 



A d 3 



q 



k° kq k 2 



q • k ' q 2 



~, ■ (8) 



To extract the divergent contribution one can 
again expand in powers of k • q/q 2 and k 2 /q 2 , 
but the dependence on the first argument of / 
remains. Therefore the linearly divergent part of 
IT 4 " has a non-trivial dependence on fc°/q ■ k, 



9 Z TA / / 

Jq 



k° 



, 0, + UV finite. (9) 



(I've been somewhat sketchy, for the precise form 
of the divergent contribution on the lattice see 
fll] Jf^] . ) Such a divergence is clearly a disaster 
because it is non-local in space and time.[] Effects 
of these divergences in lattice simulations have 
first been observed in Refs. [[44] , ff5|| . 

It was argued in jl6) that the UV divergent 
part of IP" has the same qualitative effect as the 
hard modes in the quantum theory (see Sec. 4.2), 
which is that it slows down the dynamics of the 
soft gauge field modes. In the continuum limit 
A — > oo the characteristic time scale of the soft 
modes should diverge, so that the Chern-Simons 
diffusion rate Tcs should go to zero. This be- 
haviour was indeed observed in a recent calcula- 
tion on a very fine lattice [^7j. 

4.2. Effective classical theories 

If the high momentum modes, which do not 
have large occupation numbers and are thus not 
classical, are so important for the soft field dy- 
namics, can one still make use of the classical 
field approximation? This is indeed possible. The 
modes with momenta larger than g 2 T are weakly 
coupled and can be integrated out using pertur- 
bation theory. One finds that the correct effec- 
tive classical theory, at leading order in g and 
log(l/g l ), is described by the equation of motion 

m 

DxB = 7 E + C (10) 

instead of Eq. (||). The constant 7 is propor- 
tional to T/log(l/g). Furthermore, £ is a Gaus- 

6 In QFT theory the linear divergences in Eqs. (|?j), ^ are 
cut off by the Bose-Einstein distribution function. 



sian white noiseQ Its expectation value vanishes 
and it is entirely determined by its 2-point func- 
tion 



{Q{x)d{x'))=2T 1 8 ab 5 l] 5^{x-x'). 



(11) 



Here a, b are indices of the adjoint representation 
of the gauge group. Recently the leading log ef- 
fective theory ( [To| ) was extended to include the 
Higgs field @. 

The physics behind Eq. (|l0|) is easily under- 
stood. The soft gauge fields evolve in time, which 
means that the color electric field is non-zero. 
The hot plasma is a conductor and the electric 
field induces a current which is carried by charged 
particles corresponding to the hard (k ~ T) field 
modes. The particles suffer collisions with a typ- 
ical mean free path of order (\og(l / g)g 2 T)~ 1 . In 
the weak coupling limit log(l/<?) ^> 1 this is small 
compared to the wavelengths of the soft non- 
perturbative modes. Thus the soft modes can- 
not resolve the trajectories of the particles and 
the induced current is simply a number times the 
electric field. The noise term is due to thermal 
fluctuations of the field modes which have been 
integrated out. 

Eq. ( ]i"o| ) is similar to the equation of motion of 
an ordinary (abelian) magnetic field in a conduct- 
ing plasma, where the role of 7 is played by the 
electric conductivity. (In analogy with electrody- 
namics one calls the coefficient in Eq. ( |i"o| ) color 
conductivity) . The key difference is that the elec- 
trodynamics analogue of Eq. (|l0|) is valid to arbi- 
trary accuracy if one considers fields with longer 
and longer wavelengths. This is possible because 
abelian magnetic fields are not screened. On 
the other hand, non-abelian magnetic fields are 
screened on a length scale of order (<7 2 T) _1 . Thus 
there is nothing like a magneto-hydrodynamic 
limit in hot non-abelian plasmas. In fact, Eq. ( [Toj ) 
is only valid at leading and, as was shown re- 
cently in Ref. [|lj, at next-to- leading order in 

loga/s)- 1 . 

From Eq. ([l0]) one can estimate the characteris- 
tic frequency ko of the soft modes as follows. The 
rhs is of order (<? 2 T) 2 A, since it contains two co- 
variant derivatives of the vector potential A. The 



7 "White" refers to the frequency spectrum of £ being flat. 
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electric field which contains one time derivative 
can be estimated as fcoA. Thus the characteristic 
frequency is 



fco ~ \og(l/ g)g i T. 



(12) 



This is much smaller than the estimate kg ~ g 2 T 
which one would obtain assuming that high mo- 
mentum modes decouple from the soft dynamics. 
From Eq. ( [l2| ) one can estimate the Chern-Simons 
diffusion rate, which is the number of topological 
transitions per unit time and unit volume, as 



r cs ~fc |k| a ~log(l/ 5 )a 5 T 4 
where a = g 2 /(in). 



(13) 



Eq. ( |12| ) was obtained through a sequence of ef- 
fective field theories by integrating out more and 
more high momentum degrees of freedom. The 
first is the so called 

Hard thermal loop effective theory 
It is the result of integrating out "hard" physics 
associated with k 2 of order T 2 , and it can be de- 
scribed by the classical field equations of motion 
53 



D E + D x B 
D E 

v-DW 



'D 

= m\W 
= v • E. 



i) 



(14.a) 
(14.b) 
(14.c) 



The field W — W(x,v) represents the fluctua- 
tions of adjoint color charge due to the hard de- 
grees of freedom, which act like particles with 3- 
velocity v, v 2 = 1. Furthermore, 



(l,v). 



(15) 



The 4-current on the rhs of Eqs. ( 14.a|) and ( 14.b| ) 
is given by 



d 2 fl 



4tt 



(16) 



times the square of the leading order Debye mass 
tod oc gT. Originally hard thermal loops where 
discovered when trying to develop a consistent 
perturbative expansion of Greens functions with 
external momenta of order gT 52 1 . The effect of 
hard thermal loops on the soft dynamics was first 
realized in Ref. (4^] where the fields in Eq. (14) 
where interpreted as being soft. Then the main 



effect of the hard particles is Landau damping |54| 
of the soft modes, resulting in the estimate fco ~ 
g 4 T, which is the same as Eq. (|l2|) but without 
the logarithm. It was then realized that the fields 
in Eq. ( |l4[ ) must be interpreted as modes with 
both k - gT and k - g 2 T g§. 

Non-abelian Boltzmann equation 
The k ~ gT modes in Eq. ( |l4| ) are weakly cou- 
pled and can be integrated out in an expansion 
in g. They are responsible for small angle scat- 
tering between the hard particles. They induce a 
collision term and noise term in Eq. (14.c), and 
one obtains [E8| 



D E + D x B 
D E 

(C + v- D)W 



v E 



W° 



(17.a) 
(17.b) 
(17.c) 



where now the gauge and W fields contain only 
spatial Fourier components of order g 2 T (modulo 
logarithms of 1/g). The linear collision term CW 
breaks time reflection invariance and describes 
dissipation caused by the modes which have been 
integrated out. These modes also perform ther- 
mal fluctuations which pump energy back into 
the g 2 T modes through the Gaussian white noise 
£. (For other approaches, which do not make 
use of the hard thermal loop effective theory, see 
Refs. H-||.) 

Langevin equation(s) 
The collision term C in Eq. ( 17. c ) is of or- 



der log(l/ g)g 2 T, while derivatives of the non- 
perturbative modes are of order g 2 T. At leading 
order in log(l/g) _1 the latter can be neglected. 
Then the current on the rhs of Eq. (17.a) is a lo- 
cal and linear function of E and £. Due to the 
strong damping one can neglect the kinetic term 
L» E in Eq. 0^ ) which leads to Eq. @. 

Starting from a non-local Langevin equation 
which is valid to all order in log(l/<?) _1 |5(]] it 
was shown that the logarithmic approximation 
can be systematically improved in an expansion 
in log(l/p)- 1 Q. It was found that Eq. @ is 
still valid at next-to-leading order in log(l/g) _1 
if one includes a next-to-leading log correction in 
the color conductivity 7, which was computed in 
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Figure 3. Chern-Simons diffusion rate 

parametrised by k' (see Eq. (Q) vs the 
magnetic screening scaie over Debye mass (for 
details see text). 



4.3. Lattice results 

The effective theory ( |io|) is perfectly suited for 
lattice simulations because it is UV finite [p5| . It 
was used to calculate the Chern-Simons diffusion 
rate in Ref. Due to Eq. ( |l3[ ) the result can 

be parametrized as 



Tcs = k 



a 2 T 2 

19 r a 5 T 4 



with 



(18) 



(19) 



The dimensionless constant n was determined as 



10.8 ±0.7. 



(20) 



Near a weakly first order electroweak phase tran- 
sition one has to take into account the Higgs field 
which requires an extension of Eq. ([To]) jl9| , the 
resulting effect on Tqs is of order 20%. 

The hard thermal loop effective theory (|lj) is a 
classical field theory which correctly incorporates 



the effect of the hard modes. Unlike the Langevin 
equation (|l(]) its validity is not restricted to the 
logarithmic approximation, but it does not have 
a continuum limit. |^| It was argued that it should 
still be useful for lattice simulations as long as the 
lattice spacing is not too small frJofl . 

Two different lattice implementations of the 
hard thermal loop effective theory have been de- 
veloped. In Ref. |6lJ the hard modes are rep- 
resented by classical charged particles. This 
method was further developed and used to com- 
pute r cs in Ref. |6(|. Starting point of Ref. || 
were the equations of motion ([l4]) . The W fields 
were expanded in spherical harmonics and only 
components Wi m with I smaller or equal to some 
cutoff l max were kept. 

The results are shown in Fig. 3. The rate Tcs is 
parametrized by k' through Eq. (18). The black 
symbols are the results of Ref. ]62| for 3 different 
values of the lattice spacing a, where 

4 



0o 



9 Ta 



(21) 



The open circles are the results of Ref. |i0| . Also 
shown are the results of Ref. [|7j for the classical 
Yang-Mills theory (|j). f\ The different methods 
give surprisingly similar results. The dashed line 

is a linear fit k' = C\ + C2 f • The full line 

is a fit which includes the leading log term, n' = 
tlog(p^) +c. 

The results shown in Fig. 3 are about a factor 
5 larger than the leading log result. However, 
including the next-to-leading log correction to the 
color conductivity 7 increases the value of k' in 
Eq. ( |l^) by about a factor 5 jsijl , which gives a 
remarkable agreement with the results in Fig. 3. 

5. SUMMARY 

At present we see an interesting interplay be- 
tween analytic and numerical lattice work on non- 

8 In th^UV the effect of the W fields can be ignored, and 
Eq. ( |14| ) has the same UV problems as classical thermal 
Yang-Mills theory. 

9 For the results of Ref. jfj the value of ; n E qs jigj ) 
and (1191 was chosen such that the classical polarization 
tensor7H), averaged over the directions of k, has the same 
fco — » asymptotics as the p hysi cal hard thermal loops, 
following a suggestion in Ref. |43l. 
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equilibrium field theory. It is motivated by the 
physics of the early universe and of heavy ion col- 
lisions. Real time simulations of non-equilibrium 
quantum field theory are not possible. A common 
tool is the classical field approximation which 
can be applied if the field modes of interest have 
large occupation numbers. There are several ap- 
proaches which go beyond the classical field ap- 
proximation. They require other approximations, 
such as Hartree-Fock, large N. It is possible to 
use the classical field approximation even when 
quantum effects are important by constructing 
effective classical theories. These can differ sub- 
stantially from the naive classical limit of the un- 
derlying quantum field theory. 
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